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FOREWORD 


This  report  was  prepared  by  the  Industrial  and 
Experiment  Station,  Department  of  Mechanical  Engi¬ 
neering,  University  of  Florida,  Gainesville,  Florida, 
under  Contract  No.  FC8635-71-C-0073  with  the  Air  Force 
Armament  Laboratory,  Eglin  Air  Force  Base,  Florida, 
during  the  period  from  9  December  1970  to  9  December 
1971.  Lieutenant  Robert  J.  Earner  (DLWG)  monitored 
the  project  for  the  Armament  Laboratory. 

The  principal  investigator  for  the  contractor 
was  Dr.  J.  Mahig. 

This  report  consists  of  two  volumes.  Volume  I 
is  devoted  to  the  Six-Degree-of-Freedom  Simulation 
while  Volume  II  is  concerned  with  the  Three-Degree- 
of-Freedom  Simulation.  This  is  Volume  II. 

This  technical  report  has  been  reviewed  and  is 


ABSTRACT 


This  report  describes  in  detail  a  three-degree-of- 
freedom  program  which  can  be  used  to  determine  the  tra¬ 
jectory  and  miss  distance  of  a  guided  missile  system. 

The  options  for  the  program  are  such  as  to  permit  varia¬ 
tion  of  the  aerodynamics,  seeker,  autopilot,  actuator, 
and  missile  motor  performance  for  the  purpose  of  accu¬ 
rately  simulating  a  given  missile  design  and  evaluating 
the  effects  of  any  changes  in  system  parameters.  Suffi¬ 
cient  detail  has  been  included  in  the  text  to  minimize 
the  effort  needed  to  update  or  modify  the  program  pre¬ 
sented. 
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SECTION  I 
INTRODUCTION 


purpose  of  the  report  is  to  describe  in  detail 
the  equations  used  and  the  form  required  for  the  input 
data  to  a  three-degree-of-f reedom  simulation  of  a  laser 
guided  missile  system.  The  primary  purpose  of  a  three- 
degree-of-f  reedom  simulation  is  to  make  possible  the 
study  of  the  characteristics  of  a  given  missile  system 
quickly  and  economically .  Since  cross-coupling  effects 
influence  a  missile's  performance  through  the  seeker, 
autopilot,  and  missile  dynamics,  any  results  obtained 
must  be  considered  as  preliminary  until  justified  by  a 
six-degree-of-f reedom  simulation . 

Examples  of  some  of  the  studies  which  can  be  use¬ 
fully  made  are  optimization  of  the  seeker,  autopilot 
gain,  and  noise  sensitivity.  The  noise  sensitivity 
studies  include  not  only  electronic  noise  but  spot  mo¬ 
tion  and  signal  loss  simulation.  Once  the  parameters 
for  a  given  missile  have  been  entered,  it  can  be  used 
efficiently  even  on  a  small  computer  since  the  program 
only  requires  17,000  words  of  core. 

The  program  for  a  specific  mi^ile  is  given  in 
Appendix  I.  This  program  is  made  up  of  four  parts: 
the  main  program  and  the  subroutines  DER,  MODAMS,  and 
GAUSS.  The  main  program,  modifies  input  data  to  conform 
to  program  logic  and  maintains  sole  control  for  batch 
processing  requirements.  The  subroutine  DER  retains 
control  only  over  a  given  missile  flight.  The  purpose 
of  this  routine  is  to  determine  the  current  value  of 
the  aerodynamic  coefficients  and  angle  of  attack;  pro¬ 
vide  limiters  (e.g.,  fin  deflection  and  electronic 
saturation) ;  and  determine  the  present  value  of  the 
derivative  of  the  state  variables.  The  subroutine  DER 
then  updates  the  state  variables  by  calling  the  subrou¬ 
tine  MODAMS  which  integrates  the  entire  system  of  equa¬ 
tions  using  an  Adams-Moulton  predictor-corrector  method 
with  a  Runge-Kutta  start.  The  subroutine  GAUSS  is  used 
by  DER  to  generate  a  3et  of  random  variables  which  are 
used  to  determine  the  position  of  a  laser  spot  and  the 
apparent  location  of  the  target  with  respect  to  the  mis¬ 
sile.  The  routine  is  also  used  to  determine  pulse  loss 
through  the  variable  VLAZRP. 
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All  the  aerodynamic  coefficients  are  either  a  func¬ 
tion  of  Mach  number  or  of  the  fin  deflection.  The  aero¬ 
dynamic  coefficients  used  in  this  program  are  plotted  in 
Appendix  II.  These  plots  are  used  to  develop  the  aero¬ 
dynamic  coefficients  as  approximate  functions  of  the  de¬ 
pendent  variable  Mach  number  and  fin  deflection.  The 
procedure  makes  DER  compute  much  faster  than  would  be 
possible  through  a  table  look-up  procedure  and  need  not 
detract  from  the  accuracy. 
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SECTION  II 


STATE  VARIABLES 

A  discussion  of  state  variable  techniques  is  given 
in  many  control  theory  books ' 1 • 2 •  3 ) .  However,  a  dis¬ 
cussion  of  these  techniques  is  appropriate  at  this  point 
since  they  are  used  in  the  program  to  determine  the  mis¬ 
sile's  flight  path. 

The  equation  of  motion  governing  pitch  of  the  mis¬ 
sile  i3  given  below  in  the  familiar  form: 

lyy0  =  M2  (1) 

where  m2  is  considered  i  function  of  time  and  the  state 
variables.  To  convert  this  to  state  variable  notation, 
the  second  order  equation  must  first  be  converted  into 
two  first  order  equations.  (This  is  also  true  for  an 
nth  order  system  which,  in  order  to  be  converted  into 
the  itate  variable  form,  must  first  be  converted  into 
n  '  order  systems.)  Equation  (1)  is  converted  through 

*  owing  definitions: 

6  =  02 

e2  =  m2 

Now,  since  M2  is  a  function  of  8  and  6  ana  t,  then  the 
system  may  be  defined  in  state  variable  form  as 

0  =  82 

8 2  -  M2  ( 6 , 8  j  ,t)  . 

This  same  process  can  be  used  for  an  nth  order  system. 

A  list  of  the  state  variables  and  their  derivatives 
in  DER  is  given  in  Table  I . 


TABLE  I.  TROGRAM  VARIABLE  AND  STATE  VARIABLE  IDENTIFICATION 


State  Variable 

State  Variable 
Derivative 

1. 

X  -  Y (1,1) 

X  -  YP (1,1) 

2. 

X  -  Y  (2 , 1) 

X  -  YP (2,1) 

3. 

Z  -  Y  (3, 1) 

Z  -  YP (3,1) 

4. 

Z  -  Y(4,l) 

Z  -  YP  (4,1) 

5. 

Xs  -  Y  ( 5 , 1) 

As  -  YP (5,1) 

6. 

X  -  Y (6 ,1) 

X  -  YP (6,1) 

7. 

<5  -  Y(7,l) 

6  -  YP  (7 , 1) 

8. 

6  -  Y (8,1) 

6  -  YP (8 , 1) 

9. 

6  -  Y  (9,1) 

6  -  YP (9,1) 

10. 

6  -  Y (10 , 1) 

6  -  YP  (10,1) 

11. 

0  -  Y (11,1) 

6  -  YP  (11,1) 

12. 

R  -  Y (12,1) 

R  -  YP (12,1) 

13. 

4>  -  Y (13,1) 

4>  -  YP (13,1) 

14. 

4>M  -  Y  (14 , 1) 

*M  -  YP (14  , 1 ) 

Note:  State  Variables  R,$,4>M  are  the  results  of  propor¬ 
tional  guidance  implementation.  The  implementation  is 
assumed  ideal . 


r 


SECTION  III 


SYSTEM  EQUATIONS 
3.1  Equations  of  Motion 

The  equations  of  motion  for  the  three-degree-of- 
freedom  simulation  used  in  the  program  are  developed 
as  follows: 

The  forces  in  body  axes  - 

F,  =  -(CA  +  2CA6|6j  )qS  +  TH 

F:  “  -(cNaa  +  cN5fi>9s 

M2  =  (-CNa(Xcp  -  Xcg)a  +  (XCg  -  Xc)6)qSd 
-CM06(pVA/4)SdJ 

Table  II  defines  the  relationship  between  the  com¬ 
monly  used  aerodynamic  coefficients/  used  above,  with 
the  variable  names  used  in  the  fortran  program  listing 
shown  in  Appendix  I. 

The  equations  of  motion  in  the  earth  fixed  coordi¬ 
nate  system  are : 

mX  =  FjCose  +  F3Sin6 

mZ  =  -FjSinS  +  F3Cos6  +  mg 

I Y  Y  ^  =  ^2 

where  the  quantities  a,  V  ,  and  q  are  defined  as  follows: 

a  ~  vza/vxa  <  150  * 

VA  =\fVxA2  +V2A2 
q  =  *SPVA2 

and  the  variables  VXA  and  VZA  (the  velocity  along  and 
perpendicular  to  the  missile  centerline)  are  given  as: 

VXA  =  *  Cos6  "  ^  Sine 

VZA  =  X  Sine  +  Z  Cose 
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TABLE  II.  EQUIVALENCE  OF  AERODYNAMIC  COEFFICIENT  NOTATION 


Aerodynamic  Coefficient 

Prograun  Variable 

CA 

CA 

% 

DCA 

CN 

a 

CNA 

% 

CND 

c 

M0 

CMT 

% 

CMD 

Figure  1  and  Table  III  respectively,  contain  the 
coordinate  system  used  in  the  derivation  of  the  equa¬ 
tions  of  motion  presented  above  and  define  those  other 
terms  used  as  necessary. 


Figure  1.  Coordinate  System 
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TABLE  III.  VARIABLE  LISTING 

X  -  Position  of  missile  in  X  direction 

Z  -  Position  of  missile  in  Z  direction 

Xs  -  Line  of  sight  of  seeker 
X  -  Line  of  sight  of  missile 
6  -  Fin  Deflection 
6  -  Pitch  of  missile 

tg  -  Guidance  delay 

m  -  Mass  of  the  missile 
Iyy  -  Moment  of  inertia 
q  -  Dynamic  pressure 
S  -  Reference  area 
d  -  Diameter  of  missile 

p  -  Density  of  sir  (slugs  per  cubic  feet) 

TH  -  Thrust  force 
g  -  Gravity 

V^  -  Velocity  along  missile  axis 

VZA  ”  Velocity  perpendicular  to  missile  axis 

V^  -  Total  velocity  of  the  missile 

a  -  Angle  of  attack  of  the  missile 

Xca  -  Initial  location  of  the  center  of  gravity 
^  (calibers)  (Equivalent  to  missile  diameters) 

X  -  Location  of  center  of  pressure  (calibers) 
cp 

Xc  -  Location  of  fin  (calibers) 
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3.2  Seeker  Equations 

The  position  of  the  target  is  always  assumed  to 
be  at  the  location  of  the  origin  of  the  coordinate 
system.  Thus,  the  line  of  sight  may  be  given  as 

X  =  tan1  (2/X)  . 

If  the  seeker  is  assumed  to  be  a  PLG  seeker,  then 
the  variation  of  the  seeker  axis  in  terms  of  the  mis¬ 
sile  orientation  may  be  given  as 

X£  =  -6(1  +  tsS)  . 

The  angle  between  the  line  of  sight  and  the  seeker 
axis  may  be  given  as 

e  =  X-Xs. 


3 .3  Autopilot  Equations 

The  block  diagram  for  the  autopilot  used  in  the 
program  is  given  below  in  Figure  2. 


Time  Sample  and 
delay  hold  _ 


..  I  +  a„TflS 

Sc 

1 

8 

Ka  l  +  T,s 

1  +  TJT  s  *  "ijl 

Kfl-  9.0  sec  ^  “  .  5 

Tg-.08»ec  a>  -  20  rad/eec 

a8  "  lO.Oeec 
t| _  .25  tee 

Figure  2.  Autopilot  Clock  Diagram 
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3 . 4  Nominal  Thrust  Profile 
The  initial  thrust  is 

7500  lb;  0  <  t  <  0.8  seconds 
which  then  linearly  decreases  to  zero  at  1.3  seconds. 

3.5  Current  Input  Configuration 

The  input  data  for  the  missile  configuration  used 
is : 


d  =  0.416 
S  =  0.136 

XCg  =  10.6  (initial)*  [8.75  at  burnout] 
Xc  =  8.4 

M  =  4.21  (initial)*  [3.48  at  burnout] 
Iyy  =  29.56  (initial)  [25.04  at  burnout] 
TS  =  5.0. 

3,6  Guidance  Laws 


Any  guidance  law  may  be  simulated  by  the  program; 
however,  for  convenience,  some  guidance  implementations 
are  prepackaged  into  the  program.  Thus,  at  the  beginning 
of  the  subroutine  DER  is  a  list  of  commented  cards  con¬ 
taining  the  following  statements: 

(a)  Pursuit  Guidance 

(b)  PLG  Guidance 

(c)  Sidewinder  Guidance  (PNG) . 

These  statements  head  the  required  changes  in  the  pro¬ 
gram.  needed  to  implement  that  form  of  guidance.  Thus,  for 
example,  implementation  of  PLG  guidance  only  requires  the 
replacement  of  the  current  calculation  of  YP(5,1)  by  the 
one  following  the  card  with  the  statement  PLG  GUIDANCE. 


•Variation  in  XCg  and  M  during  motor  burn  is  assumed  pro¬ 
portional  to  the  ratio  of  the  impulse  imparted  to  the  total 
impulse  available  from  the  motor. 
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SECTION  IV 


PROGRAM  DESCRIPTION 


4.1  Program  Aids 


In  order  to  effectively  coordiante  the  variables 
used  in  the  program  with  the  system  equations  developed 
above,  Table  I  is  to  be  used  for  the  identification  of 
the  state  variables  and  Table  IV  will  provide  defini¬ 
tions  for  the  other  variables  used  by  the  program.  In 
addition.  Figure  3  defines  the  form  in  which  the  input 
must  be  prepared  and  the  output  format  of  the  program. 

A  complete  listing  of  the  program  and  its  output 
is  given  in  Appendix  I.  Appendix  II  contains  the  aero¬ 
dynamic  curves  for  an  actual  missile.  These  aerodynamic 
coefficients  are  all  given  as  a  function  of  Mach  number. 
Since  accurate  functional  representation  of  some  of  these 
curves  are  very  complex,  simple  mathematical  expressions 
were  chosen  which  are  true  only  over  defined  Mach  number 
regimes.  The  adequacy  of  the  representation  may  be  as¬ 
sessed  by  the  reader  by  comparing  the  results  obtained 
from  the  equations  listed  in  the  program  with  the  aero¬ 
dynamic  curves  they  are  intended  to  represent.  These 
curves  are  shown  in  Figure  II-l  through  Figure  II-8. 
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Input  Format 

FORMAT  (1P6E12.4) 

FPT, H, TEN , AR, DVI , VI ,DTH, ALP,  RANGE, LOS, GBIAS, 
VLAZRP , XLASR, DSPEED ,  TZ 

Y(l,l),  Y  (  2, 1)  ,  ¥(3,1),  Y  (4,1) . Y(6,l) 

Y(7 , 1) . Y  (12,1) 

Y  (13, 1)  ,  Y  (14,1) 

Output  Format 

FORMAT  (1P9E12.4) 

Y  (1 , 1}  ,  YP  (1,1)  ,  Y  ( 2 , 1 )  ,  YP  ( 2 , 1) .  (5,1) 

YP  ( 5 , 1)  ,  Y  (6,1) . YP  (9,1) 

Y  (10, 1)  ,  YP  (10,1) . YP  (13,1)  ,  Y  (14,1) 

YP  (14,1) 

X,TACC 

Figure  3.  Input  -  Output  Format. 
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PROGRAM  NOMENCLATURE 


TABLE  IV. 


H 

TEN 

AR 

DVI 

VI 

DTH 

ALP  ! 
range! 

los  ; 

GBIAS: 
FPT  : 

i 

TIM  ! 
■ 

THF  j 

i 

MAC  : 

i 

ES 

I 

XCP  j 

j 

D 

| 

TACCi 


Time  increment,  DT 

Angular  orientation  of  thrust  vector  t 

Offset  of  rocket  motor  from  center 
line  of  missile 

Change  in  missile  velocity  in 
feet/second  (initial  condition) 

Initial  velocity  of  missile 

Increment  in  the  line  of  sight 
(for  next  run) 

Angle  of  attack  (in  degrees) 

Range  of  target  along  line  of  sight 
(feet) 

Line  of  sight,  angle  in  degrees 
(always  positive) 

Gravity  bias  term 

Number  of  time  increments  not  to  be 
exceeded 

Current  time,  initial  time 

Time  of  next  laser  pulse 

Thrust  factor  (not  implemented) 

Mach  No. 

Surface  area 

Location  of  c.g.  in  missile  diameters 

Missile  diameter 

Total  acceleration  magnitude 

Angle  of  attack  a 


AL 


TABLE  IV.  (CONCLUDED) 


IXX 

Moment  of  inertia  in  X  direction 

IYY 

Moment  of  inertia  in  Y  direction 

IZZ 

Moment  of  inertia  in  2  direction 

CA 

Drag  coefficient 

DCA 

Induced  drag  j 

CND 

Control  vane  normal  force 

CNA 

Body  normal  force 

CMD 

Tail  misalignment  coefficient 

XCP 

CP  location  -  diameters  from  nose 

CMT 

Pitch  damping  coefficient 

IC 

Total  number  of  increments  of  initial 
condition  DTH 

IK 

Printout  occurs  every  IKth  increment 
of  H 

N 

Number  of  state  variables 

XLASR 

Standard  deviation  of  laser  spot  in 

X  direction 

VLAZRP 

Percent  of  time  laser  pulse  acquired 
by  seeker/100  (used  in  pulse  loss 
logic) 

TZ 

Length  of  time  for  which  seeker  is 
caged  (pursuit  guidance) 

TACC 

Total  accleration  of  missile 

DSPEED 

Counter  (if  greater  than  one,  aircraft 
speed  is  incremented  by  DVl  for  next 
launch 

DTH 

Increment  by  which  initial  orientation 
is  changed  for  next  launch 

APPEN1*  X  I 
PROGRAM  LISTING 
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EP*0. 

CCPKN/VAP/TEN,AR,TtM,TFF,.1TK,  TACC ,GBI AS 
•,YFl ,VLA2RP,XLASR 
•,EP,FS,ER,TZ 

CIPEASICAY(U,61tYPtt*,6),VIH) 

EXTEPAALCER 

REAUQS 

1  CCMINU6 

X*0. 

T  I  KiO • 

IC*10 
I  T  *0 
i«*e 

RE  AC  (5  ►<500IFPT,H,TEN,  AR , OV 1 »V 1,DT H .AL P ,RANGE . LOS , G8 l AS 
• , VLA2RP ,XIASR 
•.CSPEEC.T2 

WPtT£<6,900]FPT,H, TEN, AR.0V1, VI, OTH»ALP, RANGE, LOS, GBIAS 
•,VIA2RP,XIASR 
•.CSPEEC,T2 
900  FC-RPAM  1P6E12.4) 

PPT  =  FPT 
C  T  HaCT  P  /  5  7. 6 
CCAI-l  ,  IA 
A  Ylt.l)*Ok 

Y (  12,1 I»l000v 
eBC=;[C/2«CTK 

2  CONTINUE 
vi=vi»cvl 
PV=LCS/57.6 

Y( 1 , 1 1  =  -R ANOE»CCSC  PV 1 

Y( 3,II*-RANCE»S1NIPVJ 

Y(6,l)*PV 

Y(5,U  =  Y(6,n 

Y  (  10,1  )*-Y«6ilMB6C 

BBC*0. 

PT=-Y( 10.1)«ALP/57.6 
Y(2,l)-Vl«CCS(PTJ 
Y(4,i)«Vl«SINIPTI 
FS*0. 

£R-0. 

EP=0. 

UP  I  TE  (6  ,920  )  (Y<K»l  >,K»1,N| 

CC5l=l,l4 

5  VUI-Yl  III) 

6  CONTINUE 

CALLPCCAHSI  X1H.PPT.N,  V.YP,  CERi  1K1 

X*0. 

TlP^O. 

IF!CSPEEC.GT.1.)G0TU2 
I  T  -  I  T  ♦  1 
CC10I*1  >14 
10  Y(I,tl>Vttl 

Y(  10.1  ,*Y(10.1I-CTMIT 
Y(5,II-YI10U) 

V(S,II*-YI10U» 

IFIT7.GT.0. ) GO  T  0 1 2 
Y(5,l»-VC6, II 
12  CCMINUE 
V2  »Yl 

Y(2,ll»V2»CCS(Y  10,11) 

Y(«,l)  =  -V2«31NCYIlC,in 
EP*0, 
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U  <->  O  <->  U  UUUUUUUUU 


FS*0. 

£R*0. 

TIFsO. 

WRI TE <6,900 IFPT.Fi.TEN,  AR , OV l . V 1 ,0TH, ALP , R ANGE , 10$ , GB I AS 
•  ,  Vl*ZRP,XUSR 
RRI TEI6.920  1IY(X,1 I.K-l.N) 

I F t IC-IT)15,15,6 
15  CONTINUE 
920  FCRFAT ( 1P6E12.4) 

20  STCF 
ENC 

SU8RCUTrN£0£«U,N,V,VP,H,L2lLVI 
C I P  ENS I CNY (  14,6),YPI14,6) 

CCPFCN/VAR/TEN,  PR,  TIM,  TPF,  TTK,  TACC,G8IA$ 

•,*FL,VLA2RP»XLASR 
•,EP.FS.ER,TZ 
REAIPAS,IXX,IYY, IZ2.NAC 
PLG  GUICANCE 
YF (5, 1  )•-< EP  +  FS I/A4 
PURSUIT  GUICANCE 
YP<  5, 1  )  *-E P 
SICE  hlNCER  GUICANCE 
A8*20. *6.28318 
l*YP(10,l))*A8»*2 
C2M0./P! 

C3*. 001/PI 
0*1. /PI 
Pl*57. 295776 

YPI5.1  l«-EPLC2*<ER-FS)*C3«FS 
YP{ 9.11*0. 

YF<8.1  )»-2. *A7*A8*Y(8,  1 ) -Afi. *2* Y I  7, 1 1  *C4*  I  YP 1 5 , 1 1 
RANGE »SCRT<Y(1,1)**2*Y< 3, ll»»2» 

TC* .25 
TF*1.3 
TZ*0. 

TR  *  .8 

RF=. 072/32. 2 

V1*(Y(2»1I»*2*YI4» 1)*»2)»*.S 
C* • 5  «  RF • V 1 *V 1 

Al  =  < Y<2,1  >»SIN< Y( iO, 11)  *VI4,ll«C0SCY<  10.1)1  I  /  ( Y 1  2 ,  l  I  *COS  ( Y  ( 1 0  U I 
1I-TI4.1HSIMYI10|II)I 
IFIAL.EC.O. IG0TC42 
BL*ABS (At  I 
AL*/L/ei»ATANI8l> 

42  CONTINUE 
C* . 4 16 
XC*8.4*C 
ES*.l36 

TTCT«7500.* .8* 3750.* I TF-TR ) 

Y<  7,1 )*Y I  7, 1 1 *GB I AS* 1.5/57. 6 
PI*  12./57.6 
P  l  *24. /57.6 
PP*ABSIYI7, II) 

IF(FP.P1)4,4>3 

3  Yl  7 ,  1  I  *  Y 1 7,  )  )/PP*Pl 

4  CONTINUE 
IF(X.GT.TR)GCT030 
TH-7500. 

TTC*7500.*X 

ttottc/ttot 

GCTC33 

30  IFtX-Tf  131, 31,32 
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32  TH-O. 

TTC  +  l. 

CCIC33 

31  Th*7500.*ITF-X  >/ITF-TRI 

TrC*(TTCT-n-»(TF-X)».5l/TT0T 

33  TTC»1.-TT0 

XCCM8.75M10. 6-8. 751*7  TO  )»Q 
MAS«(3. 48*14. 21-3.481*7  TO) 

1 XX* 1.0686*1  .0914- .0686  I *T TO) 

IYY*25.04*I29.56-2  5.04I«TT0 
1  /.2*IYY 
MAC*V1/  l  ISO. 

IF  (  MAC.GE.2.  )MAC»2  . 

I F ( PAC.CT  ..  75)GCTO 10 

CA».4 

GCTC12 

10  IF  (MAC.GT.1.2)GCT01I 
C  A  *  1  P  A  C  -  •  56  1  ■  •  6  /  «  2  6 
GCTC12 

11  CA  =  l.44-MAC*»22/2.8 

12  IF (TH.EC.O. 1G0TG13 
CA*.93*CA 

13  CtKTINUe 
IF(PAC.GT  . . 92 ) GOTO  16 

DCA*(.25«KAC**3*.6  3  )  •  (  ABS  (  Y  (  7,  1 )  •  2 . 88  )  )  M*2. 2 
GCTC17 

16  CCA=(-(  MAC-2. 8  )**3*. 0375*1  MAC-2.8  )*.12*.4)*(ABS(Y(7,1)  ) >2 . 80 > 

12.1 

17  CCMINUE 

CNC*3.9-CCS  MAC* 3.  1*.  159/2.  8  )».l 
rFlXAC.GT..75)  GOTO  18 
CNA915.8 
GCTC21 

18  1F(MC.CT.1.)G0T019 

C6A  =  15.8*9.2/.36*(M  AC-.  75 ) 

CCTC21 

19  IF  1PAC.GT.1.5)GCT020 
CN  A*21 . 3 

GCTC21 

20  C7>A*(MAC>2.8)**2*. 71*20.1 

21  IF  (VAC.GT.1.1GOT023 
CMC*(MAC)**315*10.3  +  4.3 
GCTC24 

23  CMC*EXP(-MAC*1.  )*14.6*IMAC— l.)*4./1.8 

24  XC  P*l-<MAC-.6>  **2*  .6  32*14.0  2*<MAC-  .6)**4*.102)».416 
IF  (  PAC.CT  .1.4  )  GCT026 

CM  T*KAC«  *2*850.  *3700. 

GCTC28 

26  CMT*EXP(-MAC*1 .4)  *56C0.*(  MAC-1.4  (•  2150. 

28  CCNMNUE 

F  1*-1CA*2.«LCA*ABS  (  Y<  7,  l)  )  )  *Q*e$*TH 
F3*-1CNA*AL*CNC*Y(  7#  1)  l«Q*E3*:TP*SIN(TEN) 

Al-10. 

A2a  .08 
A3* 10. 

A4«4. 

A5-2./57.6 

-A6»TC 

A7..5 

A8-80. 

AS  *20. 

A9*l. 
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81*-(  l.*2.»A7**8«A2)/A2 

82*-)A8**2*A2*2.9A7*A8)/A2 

ej*-A8**2/A2 

e<,«-A8*  «2«A3«A9 

65»e*/A2 

86*  C9 • A1 

VP(  1,  1  )*Y(2.1I 

YP)3,1)9Y)9,1) 

YPI  7,  l |9Y)8 ,11 
ypis.Ust  I9,i) 

VP  I  10. 1  I •  Y I  11,11 
22*22-. I 
22*T  IMV.l-X 
IF1Z2.GT.0.  1GCH07 
TIMX 

2K 3 *  t  1 1. *200. /RANG S  )*6.  ) 

I F  I ZK3  . GT •  1 5 . 1GCT07 
£P  =  V  (  10. L) 

SHxO. 

AAA9O. 

CALL  GAUSS)  S2T, AAA, AAV, X) 

IF<YFL.CT.VLA2RPIGCT07 

SSTaC. 

SSx*XLASR«S l  ft)  Y  ( 6  » 1  I ) /RANGE 
CALL  GAUSS(SSX,SSr,SSXV,Xl 
FS  =  -SSXV/RAftG£*Y(  5,  II 
ER  * Y ( 6 , 1 1 
7  CCSTINGE 

IF(TZ  —  X)80, 01,81 
81  FS»-EP 

YP)5,1)*0. 

CCTC85 

80  YP(5, l |9-<CP*FS)/A9 

85  CCMINUE 

YP(6,  l)*l./IYCUU«.2*Y(3,n»»2)*«Ya,ll»VU,l)-YC3,l) 

1  *Y  (  2  ,  l  )  ) 

YP(9, 1  I  90  L • Y (9, l) ♦  82  *Y  (  8, 1 I  *8  3  •  Y  (  7 , l)*B5MER— FS I  *B6* l  YP 15 ,1 ) ) 

YPl  11  ,  1  l*-CPT*Y)  11.  I  )•(  RF*Vl/9.  1»ES»0»»2/IYY*Q»E  S»0/1YY»I  -CNA»  (  XCP 
1-XCG  I  •AL*CNC«(XCG-XC)*Y(7,  II  )-TP«AR,  IYY 
IF(A6S(YP(9,1I|,GT.  1CC0. 1YPC9, 1USIGNI  10C0  .  ,  YP  (  9 , 1 1 1 
YPI 2, l  )9(F1 9C0SIYI 10,1 1  I *F3*SIN(YC 10,1  I  I  I /PAS 
YP(9,  1  »9l-Fl«SIM  Y  ( 10,  1 1  I  *F3*C0SIY<  10,  II I  I  /PAS* 32. 2 
A  A  = 10 . 

YPl  12, 1 1»-V1*C0SIV  1 13, 1  l-YUA,  111 

YP)  13,  1  I -IV  l  ‘SlNIYJ  13,1  I  -  Y(  19, 1  I  )/Y<  12, 1) 

YP) 19, 1 l*AA*YP( 13,  1 1 
IF(X-TCI9Q,5C.50 
90  Y P  ( 9 , 1 1 90 , 

YP(8,1>90. 

YPI7.1  I  *0  . 

YP(19,1)*0. 

YPU3,1)*0. 

YP)12,1)»0. 

GCTC70 

50  IF(X-rC-FI60,60,70 

60  Y)  13, 1  193. 19159625/2. -Y!  6,  1 1 

Yt  19, 1  193.19 159625/2. *Y)  10,11-AL 
YI12.1  l*(Y(  lill*'2*Y(3,  1)**2)**«5 
YPI12,1I»-V1.C0SCVC13,1»-YU9,  111 
YPU3,  U'V^SIWlYl  I3.1I-V)  19,1  H/Y<  12,11  . 

YP  (  19, 1  l«AA»YPI  13,  II 
70  CCMlftUE 
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TACC*SCRT(YP(2, 1)**2*YPIA, 1 ) • • 2  > 

RE  TURN 
ENC 

SU8RCUTINEMCCAmS(X,P,PPT.N,L,F,0ER.  IK) 

ClPENSICNUt  IA,6),F<  1A,6> 

CCRPCN/XAR/TEN.AR,T  [P,THF.TTK, TACC.GBIAS 
•.YFL.VIAZRP.XI.ASR 
•.EF.FS.ER.TZ 

CAT  A  PREC1,  PR  SC  2 ,  PREC3,  PREOA 
.  /S5.0.  -59. 0.  37.0,  -9.0/ 

.  .CCRRl,  C0RR2 1  C0RR3.  CORRR 

.  /9.0,  19.0,  -5. Or  1.0/ 

P  2  A  *  H  /  24.0 
PI  A  PR6C1  •  M24 
P2  a  PREC2  •  H24 
P 3  *  PREC3  •  H24 

PR  »  PRECR  •  P24 
Cl  A  C0RR1  •  H24 
C2  4  CCRR2  •  H24 
C 3  *  CCRR3  •  H24 

CA  4  CORRR  •  H2R 
CATA  ALPHA2,  ALPHAS 
.  /.4,  .45573725/ 

.  .BETA21,  BETA31,  BETA32,  BEIAR1,  8EIAR2,  BETAA3 
.  /.A,  .29697761,  .I5e75964,  . 21 S 100R 3. 050965 16 ,  3.83286476/ 

.  , 0PSGA1 ,  0MEGA2,  OPS GA 3 ,  0PEGA4 

.  /.  17476028, -.55148066,  1.2C55356,  .17118478/ 

CALLCERtX.N.Ull.ll.FU.n.P.l.l) 

IF  ( PPT  .LE.  1  )  RETURN 
A2  5  ALPHA2  •  P 
A3  =5  ALPPA3  •  H 

B21  «  8ETA21  •  Y 
B31  »  0ETA31  •  H 

B32  «  86TA32  •  h 

B41  •  8ETAR1  •  p 
eR2  *  BETAR2  •  Y 

BR 3  -  EETA43  •  F 

Cl  4  0MEGA1  •  F 
C2  4  CPEGA2  •  H 
C3  R  CPEGA3  •  H 

CR  R  OHEGAR  •  F 
LL  4  PINO  I  3,  MPT- 1 ) 

CC  10  K  »  1,U 
CC  l  1  -  1,N 

0(1.5)  *  UII.KI  4  821  «  FI  I,K) 

1  CONTINUE 

A  *  X  ♦  »2 

CALLCER(A,N.0(1,4),F(1,5),F.5,5) 

CC  2  I  r  1 • N 

0(1.5)  4  UII.KI  .  831  •  Ft  I .  K )  4.  B32  *F(I,5» 

2  CCNTINUE 

A  «  X  4-  A3 

CALLCERIA.N.U(1,5),F(  1,6), P.5.6) 

CC  3  I  4  1,N 

U  (  1,5)  *  U  (  I ,  K  )  4;  841  •  F(l.K)  4  842  •  F(I,5)  4  B43  •  F(l,6) 

3  CCNTINUE 

X  •  X  ♦  P 

CaLL0ER|X,N»U(1,5),U(1,6).F,5,6) 

CC  R  I  4  l.N 

U(l.K4l)  ,  Utl.K)  4  01  »  F(1,K)  4  02  •  F ( I . 5 ) 

.  4  03  •  F (  1,6)  4  OR  »  U ( I ,6) 

R  CONTINUE 


P0P0013O 
P0H00140 
POH00150 
POPOO 1 60 
PQM00170 


PDM00190 


P0P00210 

P0HO02RG 

POP00250 

P0P00260 

PDP00270 

PDM00280 

P0P00290 

P0M00320 


PON00370 
PDM00380 
PDP00390 
POPOOR 10 
P0POOR2O 
PDP00430 
P0M00R40 

P0H00460 

P0M00470 

P0P00R80 

P0P00A90 

PDP00510 

P0P00520 

P0POO530 

P0P00540 

PDP00560 

POH00570 

POP00580 

PDP00590 


.20 


lc«k*i 

CAIICERIX.N.UI  l.KH  I  ,F(  l.Ki  l  l.F.LO.LOl 

10  CC  N  T I NUE  ►&M00610 

If  (MFT  .IT.  51  RETURN  M0N0O62O 

K  »  5 
KT-0 
JK.  IK 

CC40LK.5iMPI 
CC  20  J  *  l.N 


UIJ.K1  X  U(J.K-l)  ♦  PI  •  F ( J i K- 1 1  ♦  P2  •  FIJ.K-2)  MDM00660 

.  ♦  P3  •  F  I  J.  K  -  3 1  ♦  P4  •  FCJ.K-4)  PDM00670 

20  CONTINUE  KOM00680 

X  =  X  ♦  H 

CAllC6R(X,N.U<l,K),FCl,M.l-.K,K) 

0C  30  J  *  l.N  MOM00700 

U(J.K)  *  UIJ.K-l)  ♦  Cl  «  F(J,K1  ♦  C2  •  FlJ.K-t)  MDM00710 

.  ♦  03  •  FU.K-21  ♦  C4  •  FIJ.K-3)  F0MOO72O 

30  CONTINUE  K0MOO73O 


CALlDERlX.N.UI l.K) , F  C 1 .  K  I .  F  .  K »  X  1 
IF (K-6 132,33133 

32  K=»6 
CCTC37 

33  0C341T-U5 
0C34JT  =■  1  JN 

Ul JT,IT)»U(JT,  I1UI 

34  F(  J  T .  IT  )*F(  JT.  I T  4 1  > 

37  I  F ( JK-6 ) 36i 36*  39 

42  JK  *  JK* IK 

36  NR  ITEI6.UK  tUIJJ'JK  I  ,FI  JJ,  JK)  ),JJ*1,NI 
NB I T  t ( 6 , 1 1 1 X , T  ACC 
IF'l  jK*IK-6)42, 42,43 

43  JK  *  JK* IK 
39  JK=*JK-1 

IF  IUJ. 5)  >40140.12 

40  CONTINUE  P0MOO75O 

12  CONTINUE 
JK  =  JK-  I K  *  1 

NR  [TE(6,ll>  (  IU(  JJ.  JK1.FIJJ,  JK)  ),JJ*1,N) 
nRITE(6,U)X',TACC 
11  FCPFAH IX.1P9E12.4) 

RETURN  KOM00760 

ENC  MDM00770 

SUBS CUT INEGAUSS  IS. AM.V.TI 

CCNFCN/VAR/TEN.  AR.T  IM  ,  TFF,  T  TK ,  T  ACC  ,  CB 1 A  S 
•.YFL.VLA2RP.XLASR 
•.EF.FS.ER.T2 

CAT*  IK  *  37  77  77  777777B) 

CAT*  I E-7777777777770CGCCCCO8) 

C  S-TFE  RECUIREO  STANCARC  DEVIATION 

C  AN- IS  THE  RECUIREO  MEAN 

C  V— VAIUEOF  CCMPUTEC  NORMALLY  DISTRIBUTED  RANDOM  NUMBER 

IFIT.GT.O.l COT  04 

2  IF(IX)7.8,8 

7  l  X  =  -  I X 

8  IF! IX. GT. 999  1CCT03 

’  IX»1X*2«1 

A»0  . 

GCTC4 

3  IX>IX/10 
GCT02 

4  CC50I-1.12 

1FI  l.GT  .  1 1 1 X*  I Y 
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J-IX*262147 

1*«J.JNC.K 

17-lX.AdO.J 

8-FLCATI I Z 1 / 1.4359 7  396 10 
YFiJa.*NC.e 
I  X*  !Y 

50  A  =  A  *YF6 

\l=  <  A-6. )•$♦/« 

fit ICRN 

ENC 


l.OOOOEtOJ 
2.5000E  *00 
2.i000t  <00 


2.0000E-02 
O.OCOOt  *C0 
0 .QCCOt  *00 


0.00006*00 
0. 25006*03 
0.0CCQE»00 


O.C00OE«OO 
9. 75006*00 


0.00006*00 

2.5000E-01 


5.0000E*02 

8.00006*01 


OOCGc-02 


I  o 

I  c  o 


!  f\#  ^  n 

O  O  Q 


O  c  o 

i  I  I 
Ul  J  M 
—  >-l  tf  ,■* 

^  •  c  -  c 
*  -  -  c 
«»  '  *»  r  :rt 


-li  . :  »i  i 

—  •  —Iff 

\%  e  i 


o  o  o  o,  O' 

t  t  ; 

JJ  U  J.l  u> 

-I  r..  -g  ",  <; 

o  a  r>  ** 

’i  r.  'V. 

i/'  0  •*> 1  ff 


'  i—  ^  i  •  " 

'  *  *  O  C 

!  —  C 


»  !  1  f 


o  o  o  o  COCO  s. 


r* 

0  -no’-! 

—O'. 


' 

1  C  O  •_!  I 

I  I  I  I 


'«  — '  —  %  — 

o  o  ■  c  ■->  a  o 

i  »  lit 


—  **  J  —  g  —  —  *. 

C  O  •  <-»  *0  »-■  '  >  -J  v. 

I  till  I 

o  -  <.*  1  **•  *  ff  .* 

O  — 


0  4  0  cok 

o  ^  o  O  •  •  • 

•  •  • !  «  •  £  O  T>  • 
»•  n  •*  o  o  i  i 


Tun  1 1  •  i  •  j  -  j  :*  >  t  t 


/j  v***t  *b  ;  i  rc 


I** 


*\«  I  ■*>  —  M 
O  to  o  o 
!  I  I 
arid  aj  u 
7-  <  n  ^  o 
r  ■*  !»  — 


Id  Ujlu  W 
O  i'T  i 


t  ~-yj 

o  o  o  ;o 

*  » » 

M  ij  UJ  J 
'j  4  -  ,  -n 
W  O  O  v 


'j  —  rj 

o  o  o  o 

»  i  ; 

U  ul  ’Jj  'll  • 

♦ 

O  T  *• 

"3  -  <0  r.  » 


^  J  o  *• 
"  f\.  C  ^ 
iX  —  ^  ry 


I  1  1 

Ul'-iJ  U  Ui 

in.n.  ;>  « 
Ifl  i  ♦  N  iT 
O  >C*  •*>  — 
**  —  £  o 


■V  -1  j-  l%* 

o  o  jo  n 


o  ^  i*  o 

T  *-  O 
s  o'o  in 


c*  -n  —  v 

o  o  olo 

i  1 1 

IJ  \‘J  -J  ,UJ 
«  «V  «“  *1 

,*  ->  ■»  .  o 

T  O  O  !o 


a  n  *»  •»; 

-i  ♦  -o  XI  . 

o  a  c  'O  j 


*1  f  ^  o 

-I  r»  ~  «r 

—  CJ  »»  £j  *T 


u  a;  'J  J  ' 
<•  i  j  < ^ ' •*  • 
S>  >  ..  'V 
—  O  71  1  — 

j)  —  y  o 


C  C*  «•  !'•  * 

u  o  n  N  .* 

i/*  C  n  «j  .tf  ; 

■  ♦  I'#  a  o  4*  ®? 


—  —  ■/»<“*  —  —  —  r> 
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td-3C<2f.*b  10  *M)W«i  *i  -  10  lt> k  0  3  »?*9M  10  3P»6%43-  20  1*0  TO  30% 20*1  10  *ll*\ 

20  00  t  3 
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Figure  II-2.  Aero  Roll  Damping  Curve 
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Figure  II-4.  Tail  Misalignment  Coefficient 
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Figure  II-C.  Control  Vane  Normal  Force 
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